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Abstract 

An alternative formula is presented for the evaluation of the zeta func- 
tion values C(2fc) without the need for Bernoulli numbers. Our formula 
is recursive, and improves the efficiency with which we can calculate 
large values of the zeta function. 



1 Introduction 

The Riemann zeta function is highly explored function in the field of 
mathematics since it lies in the confluence of many fields like number 
theory, probability, applied statistics, etc. The proof of the Riemann 
conjecture has impact in the field of security due to it being directly 
relevant in the distribution of primes. Riemann Conjecture simply stated is 
as follows 'All non-trivial zeros of the Zeta function have real part 
equal to one-half. In this paper we have looked at the already existing 
formulation[2j[lJ of the zeta function at even integers. We have aimed in 
this paper to replace the relation below containing Bernoulli numbers. 

sv ; 2(2fc)! 



2 Motivation 



The formula presented above is a simple expressions but is used in the 
evaluation of zeta function at higher arguments both for the even and the 
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odd case (already existing formulae are shown in the compilation of [6], [7], 
[8] . Although the relation seems straightforward we would like to 
highlight that computation of Bernoulli numbers for large arguments of 
zeta function might cause a computational problem. The papers [1J,[S] 
contain more efficient algorithms for evaluating Bernoulli numbers for 
higher values of the argument. Hence instead of employing better 
algorithms to calculate the Bernoulli numbers faster and accurately we 
realized that an alternative to the Bernoulli numbers in the computation of 
C(2k) would be better. Few technical aspects with regard to Bernoulli 
numbers have been highlighted below. 

3 Bernoulli Numbers 

Bernoulli numbers are employed in many functions the more prominent 
ones being the computational aspects of Fermat's Last Theorem, Riemann 
zeta function, Taylor series expansions of tangent and hyperbolic tangent 
functions, etc. The general iterative relation used to calculate the bernoulli 
numbers is presented below. 



The main problem lies in the computation of higher values of Bernoulli 
numbers using the above recursive function. Many algorithms for the 
computation of Bernoulli numbers came up after 1995 when powerful 
computers and software like Maple, Mathematica failed to compute huge 
values of Bernoulli numbers. It was reached a conclusion that computation 
of B(1000)( which has approximately 2000 digits in the decimal 
representation) would be impossible and hence many researchers started 
coming up with algorithms for faster and better approximation and 
calculation of Bernoulli numbers. Although a program was created in 2002 
for evaluation of B(4,000,000), this has come up after 6 years of research 
work. We haven exactly looked into the computational aspects in terms of 
complexity, time, precision of the formula presented as shown in history. 
With the simple idea of replacing expression for Bernoulli numbers in the 
evaluation of C{2k), we have obtained the following result. The crux of the 
result was the use of Fourier Series to get a result from which we would be 
able to compare the infinite summation series for zeta function and the 




if r=0 
if r >1 
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Fourier series representation to derive an alternate representation for ((2k) 
which yields excatly same results. 



4 Main Result 

The general fourier expansion of any function g(x) with the period T. 

2mrx ^ 2mxx 
g(x) = A a + 2_^(A n cos— — + B n sm — — ) (1) 

n=0 

We consider our as 

g k (x) = (iT - X) 2k S.t. < X < TT 



A a = - Ax-TT^dx 

n Jo 



TT 



, ) 2fc dt 

2A: 



2A; + 1 



(2) 



2 

A n 2k = — (x — n) 2k cos (nx)dx 
n Jo 

2 r . . 2 i,sin(nx) f n sin(nx) . 2k _, 

= - (x - TxY k — - J -2A; / — -(x - vr) ife x dx] 

tt n J n 



2 -[- 2 ±[{x - |5 jf (- C ° S J^) (2fe - l)(x - vr) 2fc - 2 dx]] 

2 2fc (-vr) 2 ^ 1 + 2fc-l P _ ^-2^ 

7r n n n Jo 

22k 2fe _! (2fc)(2fc-l) vr 
^ (vr) ^-ag] 

47r 2fe - 2 fe (2ife)(2fc-l) 

2 2 A n2 k-2 
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I would like to keep a k = ^ and b k = g*lgfc-ij 



4ir 2k - 2 k (2fc)(2fc-l) 

Ai,2fc — 5 5 Ai,2fc-2 

A n ,2k = a k TT 2k + &fcAi,2fc-2 (3a) 
Ai,2fc-2 = ak-lTT 2k ~ 2 + b k -lAn t 2k-i 

bkA n ,2k-2 = bkak-iTT 2k ~ 2 + b k b k _iA nt 2k-4 (3b) 
b k b k -\A n ^ k -A = ^fc^fc-i«fe-2 7r2fc " 4 + b k b k _ib k -2A n ,2k-Q (3c) 



n— 1 n— 1 



Q &fe-i^42fc-2n = Ofe-n b k -iTT 2k 2 ™ + Ai,2fc-2n-2 (3d) 

i=0 i=0 i=0 

Summing the subequations: 

fc-2 n-l fc-2 

A n ,2k = a k-n II b k ^TT 2k - 2n + JJ &fc_; A n>2 (4) 
n=0 i=0 i=0 

Using Equation(lO) from Result 1 with (j=k-2) in equation (4) 



A n , 2k = %, 2k ~ 2 + ^^, 2k ~\-2k)(2k - 1) + ... 

(-2) k ~ 1 4 
' -(fc(fc-l)...2(2ife-l)(2fc-3)...3)- 



n 2(fe— 1) n 2 7T 2 



(fc?r2fc | -2(fc)(fc - l)vr 2fc - 2 | (-2) 2 k(k - ljjk - 2){2k - l)(2fc - 3)vr 2fc - 4 | ^ 



n 2 7r 2 n 2 (n 2 ) 2 

4 ,-2, 



;i% 2((^ : ) fc " 1 (Mfc-l)(^-2)....(2)(l))((2fc-l)(2fe-3)(2fc-5) (3)(1))) 



(5) 
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From Result 2 we get the following generalization. 

_ 8k(2k - iy.Tr 2k - 2 A (_i)j-i(fc-j + i) 

n 2 vr 2 ^n 2 ^ 2 i{2k -2j + 2)V 

8k(2k - l)!vr 2fc - 2 * (-i)J-i(fc-(j-i)) 

n 2 0- 1 )7r 2 i(2fc - 2(j - 1)! Jj 1 j 

We replace all the (j-l)'s in the above summation by j's and thus 
summation changes from j = 1 — >ktoj = 0^k — 1 



8k(2k - l)!vr 2fe - 2 t4 , -1 , k 



^n,2fc = -j 2^ 



n 2 Z-,^^ (2fc_2j)! 

Ak{2k - l)\7T 2k ~ 2 (g, -1 V7 . 2(/c-j) 



n 2 j^nW {2k-2j)\ 

2(2k)\^- 2 -1 ■ 1 

n 2 ^ " 2 7T 2 (2fc-2j - 1)! W 
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Substituting (7) and (2) into (1) we get the following. We realize that 
B n ,2k=0 since gk{x) is an odd function. 



„2fc 00 nfnuy 2k-2 i i 

<* - ^ - 2F+T + E pdjT^-(-) 

With x = 

^ 2fc ^ 2(2A:)!7r 2fc - 2 ^ . -1 . , 1 



(7r) =2TTT + ^ ^ ( 



n 2 7r 2/ (2fc-2j-l)! 

n=l j=0 y J ' 



2fc + l n 2 v n 2 7r 2 ' (2fc-2j-l)! 

, , fc-i 



ir 2 ^}^' (2k-2i-l)\^}n 2> { ) 



(2fc + i)! vr 2 ^7r 2 ' (2fc-2j-l)! n , 

For notational convenience we set Yl^Lii^V = C(2j) 

fc— l 



A: 



- = — ly 1 Vf-V +1 
)! vr 2 ^ l 7T 2; f2fe - 27 - D! ^n 2 ' 



(2k + l)\ K 2 j^v 2 ' (2fe-2j-l)! n , 



fc-2 



+ (2kTT)T " -<D?> J+1 (2t 4-l)! C(2i + »» 

oo fc-2 _ , 

(9) 



fc 
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Result 1: We try to simplify the product or summation in 

j 



IK 



i=0 

by finding an equivalent expression for it. 
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bkbk-ibk-2bk-3-—bk-j — 



[ h-i = 

i=0 

(-2k)(2k - 1) (-2k + 2))(2k - 3) (-2k + 4)(2k - 5) 



n 2 n 2 n 2 



(10) 



L^2^(k(k - l)(k - 2)(k - 3)....(k - (j - 2))(k - (j - l))(k - j)) 

(2k-l){2k-3)....(2k-2j + 1) (11) 



Result 2:Let us consider the last term of the Equation [5] 



' ((-J) fe -W - l)(k - 2)....(2)(l))(2fc - l)(2k - 3)(2k - 5) (3)(1))) 



n 2 7r 2 n 2 



{(2k - l)(2k - 3)(2k - 5) (3)(1)) x 



(fc(fc-l)(ife-2)....(2)(l)) = k\ 

(2k)(2k-2) (2) _ (2k)\ 



(2k)(2k-2) (2) 2 k ~ l (k\) 



We would like to generalize the Equation [5] by finding the j th term so as 
to refer it by an index. Hence the j th term of Equation [5] is. 



^((-^y-HHk - 1)....(* - j + l))((2k - l)(2k - 3) (2k - 2j - 3)))^ k '^ 

4 (-l)j-l 2 fc(fc-j + l) (2k-l)\ ) 2k -2 j+ 2 

n 2 TT 2K n 2 ^- 1 ) (2fc-2A; + 2)! ; 
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Conclusions 

Having checked the above result for few values of C(2k) using the above 
formula and verified with pQ 

«*> = -(E(> V-2j-i)! c(2j + 2) + WTTy ] ' 2ki - 1]t 

= 1 C(2) 
k = 2 -> C(4) 
fc = 4 -> C(6) 

The above recursive formula gives as an alternative for the regular 
Bernoulli dependent formula for the Q{2k) which yields perfect results for 
all values of k and has avoided the use of Bernoulli numbers. 



= 7r 2 /6 

691 x 2 11 12 

= ; 7T 

15! 
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